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Abstract 



We reconsider the extraction of the gluon condensates (a s G }, (g f„bcG } and the MS running quark masses m c j, from different 
. . . Ai„(Q 2 ) Moments and their Ratios by including PT corrections to order a], NPT terms up to (G 4 ) and using stability criteria of the 
7 ~ ' results versus the degree n (number of 2 2 -derivative). We explicitly show that the spectral part of the lowest moment A1i(0) depends 
strongly (as expected) on its high-energy (continuum) contribution, which is minimized for A / l„>3_4(0). Using higher moments and 
^ jthe correlations of (a s G 2 ) with (g 3 f abc G 3 ) and <G 4 >, we obtain (a s G 2 ) = (7.0 + 1.3) x 1(T 2 GeV 4 and (g 3 f abc G 3 } = (8.8 + 5.5) 
GeV 2 x (a s G 2 ), while our analysis favours a modified factorisation for (G 4 ). Using the previous results, we re-determine m c (m c ) 
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and find that the commonly used Ati(0) lowest moment tends to overestimate its value compared to the ones from higher moments 
where stable values of m c (m c ) versus the variations of n and the continuum models are reached. These features can indicate that 
the quoted errors of m c j, from A4l(0) may have been underestimated. Our best results from different high-n moments and their 
ratios are: m c (m c ) — 1261(16) MeV and mjjjhb) = 4171(14) MeV, in excellent agreement with results obtained in 01 using some 
judicious choices of ratios of moments. 

Keywords: QCD spectral sum rules, gluon condensates, heavy quark masses. 
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1. Introduction 

Non-zero values of the gluon condensates have been advocated 
by SVZ HH for non-perturbative QCD. Indeed, the gluon con- 
densates play an important role in gluodynamics (low-energy 
theorems,...) and in some bag models as they are directly re- 
lated to the vacuum energy density (with standard notations): 



E = — —r-{a s G ) . 

sat 



(1) 



Moreover, the gluon condensates enter in the OPE of the 
hadronic correlators [2 ] and then are important in the analysis of 
QCD spectral sum rules (QSSR), especially, in the heavy quarks 
and in the pure Yang-Mills gluonia/glueball channels where the 
light quark loops and quark condensates 1 are absent to leading 
order &@1. The SVZ value: 



(a s G 2 ) 0.04 GeV 4 . 



(2) 



extracted (for the first time) from charmonium sum rules [2] 
has been challenged by different authors [4-6]. Though there 
are strong indications that the exact value of the gluon con- 
densate is around (or most likely 2-3 times) this value as ob- 
tained from earlier heavy quarks M„(Q 2 ) |&@1, FESR JT^l 
and exponentia l 01 ill moments, heavy quark mass-splittings 
ifrjl and e + e~ lll5l4l9l| inclusive data. Most recent determi- 
nations from r-decay [20-22] (see however |23|[) give a value 
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'The heavy quark condensate contribution can be absorbed into the gluon 
one through the relation |2]: (QQ) = -{a S G 2 ) I (\2nmQ) + ...An analogous 
relation also occurs for the mixed quark-gluon condensate 



(a s G 2 ) (0.02 ± 0.04) GeV 4 , while some particular choices of 
M„(Q 2 ) charmonium moments give (0.04 ± 0.03) GeV 4 0]. 
Lattice calculations found large range of values 025l - l27ll . All 
these results indicate that the value (a s G 2 ) is not yet well deter- 
mined and needs to be reconsidered. 

In a previous paper [1], we have extracted, for the first time 
within QSSR, the correlation between (a s G 2 ) and {g^fabcG 3 ) 
by working with higher moments known to order a 2 and up to 
(g^fabcG 3 ). We have obtained: 



g 3 fabAG 3 ) = (31 ± 13) GeV 2 <a s G 2 ) 



or, in terms of the instanton radius: 



(3) 



p c * 0.98(21)GeV" 1 (4) 
if one uses the dilute gas instanton (DGI) model relation^: 



ffahcG 3 ) = 4 12^ 

5 p 2 



{a s G 2 ) 



(5) 



One may interpret the previous value of (g 3 f a bcG 3 ) as the one 
of an effective condensate which can absorb into it all higher di- 
mensions condensates not accounted for when the OPE is trun- 
cated at the D = 6-dimension. 

In the present paper, we shall study the effects of the D — 8 
condensates on the previous results considering the fact that 
these effects can be sizeable when working with higher mo- 
ments 10, IH IIH]. In the same time, we shall reconsider the 
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determination of (a s G 2 ) and m c j, from different M„{Q 2 ) mo- 
ments and their ratios by including corrections to order a] and 
non-perturbative terms up to (G 4 ). We shall also focus on the 
extraction of m c j, from the widely used M n =\{Q 2 = 0) mo- 
ments. 

2. Moment sum rules, stability criteria and optimal results 

Here, we shall be concerned with the two-point correlator of a 
heavy quark Q = c,b: 



i j d 4 x e-^iOlTtyx) (j v Q (0)f |0> . 



(6) 



where : Jq - Qy^Q is the heavy quark neutral vector current. 



Different forms of QSSR exist in the literature [4-6]. In a pre- 
vious [ 1 ] and in the present paper, we work with the moments|3 



- f 

J An 



R(t, ml) 

(t + e 2 ) n+l ' 



and with their ratios: 

AUG 2 ) 



n,/n+i{Q) 



M n+1 (Q 2 ) 



r n /n+2(Q ) 



M„(Q 2 ) 
M n+2 (Q 2 ) 



(7) 



(8) 



where the experimental sides are more precise than the absolute 
moments Ai„(Q 2 ). 

In the following, we shall use stability criteria, i.e. a minimum 
dependence of the results on the variation of the finite number 
of derivatives n. In practice, this minimum sensitivity is sig- 
naled by the presence of a a plateau or a minimum. 
We shall study later the effect of the QCD continuum models 
on the results. 

We shall denote by optimal result the one obtained within the 
previous stability criteria and which is less affected by the dif- 
ferent forms of the continuum models. 

3. QCD expressions of the sum rules 

The QCD expressions of the moments can be derived from the 
ones of the vector spectral function R. 
To lowest order, it reads : 



Rlo = " v z ) 



(9) 



where v = ^1 - 4m 2 g/t is the quark velocity. 

The a s correction is known exactly to 0(a s ) MM and an inter- 
polating formula has been proposed in I132.ll . 
To order a 2 , we shall use the approximate formula given in [33] 
and derived from the exact expression in ||34|36|]. 
To order a], the three lowest Mi(0) [37] and At2,3(0) mo- 
ments 1I38I1 are known analytically . Semi-analytic expressions 



of higher moments Ai„(Q) using Pade approximants [39] and 
Mellin-Barnes transform 114011 are also available. 
The gluon condensate (a s G 2 ) contribution to the two-point cor- 
relator is known to lowest order [2] and to order a s 114 ill . 
The dimension-six condensates ((g^f a b c G 3 ) and g A {uu) 2 ) con- 
tributions have been obtained by [8]. Convenient expressions 
for numerical analysis of different M„(Q 2 ) moments including 
the (g 3 f a bcG 3 ) term are given by [24]. We have checked some 
but not all of them. 

The (G 4 ) condensate contributions have been calculated by [28] 
to lowest order. The expressions of M„{Q 2 ) have been given by 
H and fl. 

In the following discussions, we shall not transcript all these 
previous long and tedious formulae which interested readers 
can found in the original papers. 



4. Experimental parametrization of the sum rules 

In a narrow width approximation (NWA) and for Q = c0: 



R(t) = 



47rlmll(f 

N 



ie) 



Q2 a 2 /_i 



J^M^r^-d^t - m 2 ) 



(10) 



where N - 3 is the colour number; and F^^ e + e - are the 
mass and leptonic width of the J /iff mesons; Q c = 2/3 is the 
charm electric charge in units of e; a = 1/133.6 is the running 
electromagnetic coupling evaluated at M^. We shall use the ex- 
perimental values of the J /if/ parameters compiled in Table Q] 



Table 1 : Masses and electronic widths of the Jjtp family from PDG10 1 42]. 



Name 


Mass [MeV] 


Tj/^ e * e - [keV] 


J/WIS) 


3096.916(11) 


5.55(14) 


<H2S) 


3686.093(34) 


2.33(7) 


i/<3770) 


3775.2(1.7) 


0.259(16) 


i^(4040) 


4039(1) 


0.86(7) 


iA(4160) 


4153(3) 


0.83(7) 


i/<4415) 


4421(4) 


0.58(7) 



We shall parametrize the contributions from V^c ^ (4.6 ±0.1) 
GeV using either: 

Model 1: The approximate PT QCD expression of the spectral 
function to order a 2 up to order (m 2 ltf given in [33] and the a] 
contribution from non-singlet contribution up to order (m 2 /f) 2 
given in 04311 . 

Model 2: The asymptotic PT expression of the spectral func- 
tion known to order a] where the quark mass corrections are 
neg lectedQ. 

Model 3: Fits of different data above the if/ilS) mass: the most 
recent fit is done in B43I1 where a comparison of results from 
different fitting procedures can be found. 



-We shall use the same normalization as 124 



3 A missprint of factor tt is in but does not affect the results. 
4 Original papers are given in Refs. 317 to 321 of the book in Ref. @|. 
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5. Test of the continuum model-dependence of the moments 

In this section, we test the model-dependence of the exper- 
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Figure 1: Behaviour of moments M„(Q 2 ) in units of (Am 2 )" X 10" +1 versus n 
for different models of the continuum as defined in section|4] a) M„(0) : Model 
1: green (continuous), Model 2 : red (dot-dashed), Model 3: bleu (dot); b) the 
same as Fig^i) but for M„(4m 2 ) for Models 1 and 2. 



imental side of the moments using the previous models for 
parametrizing the continuum (high-energy) contribution to the 
spectral function. The analysis is shown in Fig. QJ for the mo- 
ments Ai„(0) using Models 1, 2 and 3 for different values of n 
and in Fig QJ) for the moments M„(4m^) using Models 1 and 
2. One can deduce that this model dependence can be avoided 
when working with values of n > 3,4. One can also notice 
that for Ali (0), the continuum (high-energy) contribution to the 
moments is about (40-50)% of the total contribution, which in- 
dicates that the resulting value of m c from the low moments 
n < (2-3) will depend strongly on the appreciation of this 
high-energy behaviour which is not measured accurately as also 
emphasized in [43]. 



6. QCD inputs and higher gluon condensates 



From the different expressions of the PT series given in IB41I . we 
observe that, unlike Ai„(Q 2 = 0) where the coefficients increase 
approximately like n for large n (the same feature occurs for 



the arj term given in [38-40]), the ones of M n (Q + 0) remains 
(within a factor 2) almost constant though change sign from low 
to high moments. Therefore, we estimate the coefficient of the 
0(a\) term of the moments Ai„(Q 2 + 0) to be about: 



c 3|q 2 #o,« — ± c 3Iq 2 =o,«=i — ±5.6 . 



(ID 



which is larger than the estimate used in [1], where it has been 
assumed that the ratio of the a 2 over the a 3 s coefficients are ap- 
proximatively the same for each moment. 



We shall use the input values UlQiil 231: 
m c (m c ) = 1261(18) MeV, 
a s(ni c )\ nf= n - 0.408(14) from t - decays , 
a,(\jnlr) 2 = 4.5 x 10~ 4 GeV 6 from e + e~ . 



(12) 



The error in the value of a, is the distance between its value and 



the world average 11421 14411 . 



The QCD expressions of the moments are tabulated in [24] 
for the fixed order PT series up to order a 2 and including the 
(g 1 fabcG 3 ) condensate. The contribution of the (^(ifnfr) 2 D=6 
condensate is numerically negligible and has been omitted 
The contribution of the D=8 condensates can be found in |2 
and [7]. In general, one can form eight operators for the D=8 
gluon condensates: 



01 = 


= <Tr G 2 Tr G 2 ) , 


02 = 


= <Tr G VfI G pp Tr G YT Gf T ) 


03 = 


-- (Tr G Vfl G Tp Tr G v/J G Tp } 


4 = 


-- (Tr G w G Tp Tr G V T G^ P ) , 


05 = 


-- (Tr G Vfl G tlp G pT G TV ) , 


06 = 


= (Tr G vli G v »G Tp G Tp ) , 


07 = 


= (Tr G w G vp G^Gn , 


08 = 


= (Tr G Vil G pT G v ^G pT ) . 



(13) 

Using the symmetry properties of the colour indices and an ex- 
plicit evaluation of the trace, one can show that one has only six 
independent operators and the relation for N — 3 colours [45]: 



<9 5 + 20 7 = 



08 + 20 fi 



02 + -04 , 

03 + -01 ■ 



(14) 



Normalized to (G 2 ) 2 , the use of the vacuum saturation in the 
large TV-limit gives: 



Oi = 4-( 1 + 5^)' 

V6 T 6 N 2 -l ) ' 



03 



1/111 

02 ^A\A + 3W~l 

1/1 1 1 
° 4 = 4 12 + 2/v^I 



° 5 n*-i)' ° 6 4N\6 6N 2 -l) 

1 



0v - w(?-*7ih)> °*-4N\3 



TV 2 - 1 



(15) 



which indicates that only the first four operators are leading in 
1 IN, and the previous constraints in Eq. (TBI are not satisfied 
for large N. Moreover, the 1 /TV 2 corrections to these leading- 
term are also large for N — 3 in the case of O3 and O4, and 
raise some doubts on the validity of the 1 /TV-approximation. 
Therefore, a modified factorization has been proposed in ll45ll . 
where the D=8 gluon condensates have been expressed in terms 
of 02 which is not constrained. Normalized to (G 2 ) 2 , one has: 

0! = 30 6 = \ 3 = 20 4 =~ + 202 



05 = 07 



J- + ±07 

192 7*^2 



08 



48 



16 

+ 20 2 . 



(16) 



Ref. [45] estimates O2 using either its large N or its factoriza- 
tion value. Noting that the dominant contribution to the sum 
rule is due to G5, Ref. [145] notices that the factorization pro- 
posed in [28] overestimates the D=8 gluon condensate contri- 
butions. 

For definiteness, we use the following notations and values: 



p c = instanton radius introduced in Eq.© 
fac = Is factorisation of (G 4 ) , 
fac =i 0.5 = modified factorisation of (G 4 ) . 



(17) 



respectively from lEsIl and ll45ll . We also use the value of the 
scale M 2 « 0.3 GeV 2 estimated in [28], which characterizes the 
average virtual momentum of the vacuum gluons and quarks 
and which relates, using factorization, some of the D = 8 to the 
D = 6 condensates: 



5 f abc G a MV j»f v ) = -ig 4 m 2 M\ 



(g 1 f abc G a MV G^ A D a D a G c Afi ) = (g'f^G^G^MMm 
where 7° = ijry^iA" /2)if/ and D a the covariant derivative. 



7. Hunting (a s G 2 ) from higher moments M„(Q 2 ) 
As mentioned in the introduction, the gluon condensate plays a 
key role in QCD gluon dynamics like is the quark condensate 
(ifnfr) for chiral symmetry breaking. We have also mentioned 
that there is a spread of predictions of its value in the literature. 
The extraction of (a s G 2 ) in this paper is closed to the one using 
charmonium sum rules in the early literatures which follow the 
pioneer work of SVZ fl. 

In our analysis, we shall work with higher moments which are 
more sensitive to (a s G 2 ) but limit ourselves to the ones where 
the higher dimension-six and -eight condensates remain still 
small corrections such that the OPE remains valid. This com- 
promise choice eliminates the use of higher Q 2 = moments 
where their convergence has been the subject of hot debate in 
the past flll^ll^]. Instead the Q 2 + moments converge faster 
J3] which allow to work with higher n-values. In the following, 
we shall work with M n (Q 2 = 0) for n < 5, M„(Q 2 = 4m 2 ) for 
n < 11 - 12 and with M„(Q 2 = 8m 2 ) for n < 20 where the OPE 
still makes sense when using the values of the vacuum conden- 
sates given in the literature [|4j]. 

We extract (a s G 2 ) using its correlations with the D = 6 and 8 
condensates introduced above. We allow the instanton radius 
p c which correlates (a s G 2 ) and (g 3 f a b c G 3 ) to move from 1 to 
5 GeV -1 where the latter would be the value given by a dilute 
gas instanton model estimate |2]. We shall also use the relation 
of (a s G 2 ) and (g 3 fabcG 3 ) with the D=8 condensates if one as- 



sumes a factorization hypothesis rf28ll or its modified form 1451. 
Notice that, unlike [24], we fix m c , which is, at present, known 
with good accuracy, in order to give stronger constraints on the 
value of (a s G 2 ). We show the results as function of the number 
n of derivatives for Q = 4m c and 8m c and for different values 
of the QCD input parameters. 

One can notice from the Fig[2]that the effect of a] is relatively 




9 10 11 12 




12 14 



18 20 



Figure 2: Behaviour of (a s G 2 ) in units of GeV 4 versus n from M„(Q 2 ) mo- 
ments for different values of the QCD inputs. We use the central value of a s 
given in Eq. (12) . The region between the same curves correspond to the val- 
ues of C3 from Eq. lilt : a) M„(4m 2 ) : p c = 1 GeV - ', fac=l: green (contin- 
uous); p c = 5 GeV -1 , fac=l: blue (dot); p c = 1 GeV -1 , fac=0.5: red (dot- 
dashed); p c = 5 GeV -1 , fac=0.5: black (dashed); b) The same as FigUJi) but 
for M„(8m?) 



small. Much more stable values of (a s G 2 ) correspond to the 
case of a modified factorisation of (G 4 ) which sounds better 
founded from the analysis of [45] based on the l/N approach. 
Taking into account these remarks, we deduce in units of GeV 4 : 

{a s G 2 } = (4.8 - 9.2) x 10" 2 from M„ (4m 2 ), 

(5.6- 8.3) x 10~ 2 from W(„(8to 2 ) , (19) 

where we have used the Mathematica Package FindRoots, 
which we shall also use later on for deriving all the results in 
this paper. 

We consider as a final result the most precise determination 
from vVf„(8m 2 ) which can be written as: 

(a s G 2 ) = (7.0 + 1.3) x 10~ 2 GeV 4 . (20) 

This result goes in line with different claims 
14 , nil IH] that the SVZ value given in Eq. (O understimates 
the value of the gluon condensate \ This result agrees quite 
well with the one derived from the charmonium and bottomium 
mass-splittings using double ratio of sum rules (DRSR) 01411 : 



{a s G 2 } = (7.5 ± 2.5) x 10~ 2 GeV 4 
and from r-like sum rule for e + e - — > / - 



(21) 



(a s G 2 ) 



6.1(0.7) x 10" 2 GeV 



1 hadrons data [19]: 
(22) 



5 A compilation of different determinations can be found in Table 2 of fl4ll 
and in the book [4] (reprinted papers in Chapters 51 and 52). 
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Our result is more precise than the one in 112411 . using some par- 
ticular choice of moments, as, here, we have fixed the value of 
m c while in Il24ll a two-parameter fit (m c , (a s G 2 )) has been per- 
formed. Indeed, using as input the value of m c (m c ) in Eq. ( fT2l . 
one would deduce from the different figures given in [24]: 



(a s G 2 ) (3.5 - 7.5) x 1(T 2 GeV 4 



(23) 



obtained to order a 2 and without the inclusion of (G 4 ). This 
range of values is in agreement with the one in Eq. (|T9]» but less 
precise. 

8. Re-extraction of (g 3 f a b c G 3 ) and factorisation test of <G 4 ) 



Figure 3: Behaviour of the instanton radius p c in units of GeV versus n from 
M n (%m 2 c ) moments for the central values of m c (m c ) and a s given in Eq. I12t . 
The curves correspond to fac=0.5 for the factorisation of (G 4 > and (a s G 2 ) = 
0.070 GeV 4 . The region between the same curves correspond to the values of 
C'3 from Eq. II U : dashed (red) with (G 4 ) and green (continuous) without (G 4 ). 



Using the previous new informations, we re-extract the value of 
(g 3 .fabcG 3 ) firstly obtained in [ 1] using the sum rules approach. 
As remarked in [ 1 ], the moment Ai„(^m 2 ) can provide the most 
accurate value of {g 3 f a bcG 3 ). We plot in Fig [3] the value of the 
instanton radius p c defined in Eq. (0 versus the number of mo- 
ments for given values of m c , a s , (a s G 2 ) and the factorisation 
factor fac of the (G 4 ) condensates, . We have only shown in 
Fig [3] the curves for fac=0.5 because the one for fac =1 gives 
unrealistic values of p c . This can be an indirect indication that 
the value fac= 1 is less favoured, a result which supports the 1 /N 
result in 114511 . A similar feature is also signaled when extracting 
(a s G 2 ) because for fac=0.5, larger stabilities versus the change 
of n (see Fig [2]) are obtained. At the minimas of the curves in 
Fig |3] we deduce the optimal value of p c in GeV 1 when the 
effect of (G 4 ) is included: 



p c = 1.84 + 0.24„ s + 0.33„3 + 0.27 C 2 , 
which after adding the errors quadratically gives: 
p c = (1.84 + 0.49) Ge^ 1 

ig'fabcG 3 ) 



(24) 



(a s G 2 ) 



(8.8 + 4.7) GeV z 



(25) 



We consider this result as improvement of the previous result 
quoted in Eq. ©, which has been affected by the presence 
of (G 4 ) in the OPE (see the two continuous (green) curves in 



Fig [3] when (G 4 ) is not included|§) as (a priori) expected. This 

value of (g 3 fabcG 3 ) is in the range of lattice calculations in pure 

I — I 

SU(2) Yang-Mills [26], though an eventual future result for 
S U(3) is desirable. 

9. Tests of the convergence of the OPE 

We show some behaviour of the OPE using the set of param- 
eters obtained previously, namely the values of (a s G 2 ) and 
(g^fabcG 3 ) in Eqs. d20b and ( f25l ) and the one of a s in Eq. ( fT2l . 
The PT series include the coefficient C3 = -5.64 of a 3 , while 
the D — 4 condensate includes term to order a s . Representative 
expressions correspond to the moments where the optimal val- 
ues of (a s G 2 ) and (g 3 f a bcG 3 ) from Fig s|2] and [3] are obtained. 
Normalized to {4m 2 ) n x 10 4 , the M&$(4m 2 ) moments read: 



M 9 (4m 2 ) 



M w (4m 2 c ) 



1.1314 1.111 



0.4903 1.013 - 



0407 


0.090 


0.085 


m 4 


+ — T- 


+ 0- 






0.472 


0.152 


0.144 


m 4 




+ r- 


m 6 c 





,(26) 



while the Mn^eiRm 2 ) moments normalized to (4m 2 )" x 10 9 
read: 



Mu(8m*) 



M 16 (Sm 2 c ) 



2.3181 1.005 



0.7077 0.935 - 



0.503 


0.112 


0.189 


m 4 


+ — t- 


+ — r~ 


m° 




0.549 


0.159 


0.263 


in* 


+ T~ 


+ 8~ 


m° 





,(27) 



where one can see that the NP contributions become sizeable 
(the (a s G 2 ) contribution is 16-22% of the LO contribution) but 
the OPE continues to converge (the (G 4 ) contribution is less 
than 4%). Reciprocally, the relative large NP contributions have 
permitted the extraction of their size from the moments. 
We also show the PT expressions of the moments normalized 
to (4ml)" x 10 4 at fixed order: 

Mg T (4m 2 c ) = 1.1314(l + 0.601a, + 2.7a; - 5.64c^) , 

M P l o(4m 2 c ) = 0.4903 (l + 0.045a., + 1.136a; - 5.64aJ) ,(28) 

and (normalized to (4m 2 )" x 10 9 ): 

M^(Sm 2 c ) = 2.3181 (l+ 0.031a, + 0.77aJ - 5.64a^) , 

Mf 6 r (8m 2 ) = 0.7077 (l- 0.364a, -0.33aJ-5.64a^) , (29) 

where a s = ajn. One can note that radiative corrections to 
these higher moments are less than 1 1 % while it is about 30% 
in the case of Mq(Q) in Eq. OTI ) which makes the latter sen- 
sitive to the way how the PT series is organized (fixed order, 
contour improved,...) as mentioned in 04311 . 
The D=4 condensate contribution including the a s corrections 
normalized to the LO PT moments and without the overall fac- 
tor {a s G 2 )l(4m 2 c ) 2 read: 



M$ =4 (4m 2 c ) 
M?o = W) 



-329.4(1 - 0.862a,) , 
-433(1 - 1.673a,) , 
-413.8(1 -0.986a,) , 
-491.3(1 - 1.527a,) . 



(30) 



Again here, the a s corrections are relatively small which is not 
the case of Ai„(Q) as one can see in Eq. d33l . 



6 These values agree with the one obtained in QJ using some judicious 
choice of the ratios of moments ("13/14 and ("14/15. 
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10. Determinations of m c j from low moments M„<s(0) 
Low moments are widely used in the literature for extracting 
m c b where it has been argued that its QCD expression is under 
a good control due to the negligible contributions of NPT terms. 
Though this is absolutely true for vVti(O), the neglect of the NPT 
terms becomes questionable for other moments because they 
increase in the OPE as shown explicitly in Eq. (|3T1 i. The five 
lowest moments normalized to (4m 2 .)" read: 



Mi(0) 
M 2 (0) 
M 3 (0) 

Mm 

M 5 (0) 



0.8000 



= 0.3429 



0.2023 



0.1385 



0.1023 



1.300- 



1.350- 



1.287 • 



1.158- 



0.996 - 



0.0222 




0.0005 




0.001 


m 4 


+ 




+ 










ml 


0.0862 




0.0076 




0.007 


m 4 


+ 




+ 


m\ 


0.1780 




0.0368 




0.027 


m 4 


+ 


mf 


+ 


ml 


0.2815 




0.1172 




0.077 




+ 




+ 




m 4 


mf 


ml 


0.3620 




0.2959 




0.184 


m 4 


+ 




+ 


ml 



(31) 



which indicate that already for n > 2, one cannot neglect 
the non-perturbative contributions which are larger than 3 .4% 
(compared to a] > 1.7%) in the determination of m c . 
Another inconvenience of A1i(0) is the large contribution (> 
40% effect) of the less accurate high-energy part of the spectral 
function which implies a model-dependent continuum contri- 
bution or a dependence on the way the non accurate data are 
handled as discussed explicitly in Section [5] and in 14311 . 
Low Q 2 — moments are also affected by large radiative cor- 
rections which one can observe from their QCD expressions 
given in the literature 11241 I37M40I1 . To order a 3 s , the PT series 



normalized to (4m c )" read in our normalization: 

Mf r (0) = 0.80(1 + 2.39a, + 2.38a; - 5.64a'), 

Mf(0) = 0.3429(1 + 2.43a, + 6.1 la]- 7.64a'), 

Mf (0) = 0.2032(1 + 1.92a, + 6.12a;; - 10.48a^), 

Aif(0) = 0.1385(1 + 1.10a s + 4.40a] - 18.13^), 

Mf(0) = 0.1023(1 + 0.08a, + 2.16aJ-27.4a^), 



(32) 



where one can notice that the coefficient of a] grows with the 
order n of the moments, but the coefficient of a s decreases. 
The D = 4 contribution including the a s corrections normalized 
to the lowest order PT moments and without the overall factor 
{a s G 2 )l{4m 2 c ) 2 read: 

Mf= 4 (0) = -15.04(1 +2.480,), 

A^= 4 (0) = -58.49(1 + 1.05a,), 

Mf= 4 (0) = -143.6(1 -0.48a,), 

Mt*(0) = -283.4(1 -2.11a,), 

A(»= 4 (0) = -491.3(1 -3.80a,), (33) 

where, one should note that one cannot go beyond n — 5 be- 
cause the a s correction to the D — 4 contribution is larger than 
49% indicating the divergence of the QCD series as also em- 
phasized by 112411 . 

Then, we limit ourselves to use the relatively low moments 
vV(„<5(0) for extracting the running mass m c (m c ) within fixed 
order PT series and for a given set of NP parameters determined 
in the previous section. We show the results from the moments 
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Figure 4: a) Behaviour of m c (m c ) in units of MeV versus n from A1„(0) mo- 
ments for the central value of (a s G 2 ) in Eq. {20) and of a s given in Eq. ill) . 
We use fac=0.5 for the factorisation of (G 4 >. The curves correspond to dif- 
ferent models for the continuum defined in Section [5] : Model 1 : green (con- 
tinued); Model 2: red (dot-dashed); Model 3 (data fit) : blue (dot), b) Be- 
haviour of m c (m c ) versus different ratio of moments A1„(0) moments in units 
of MeV. The inputs and legends are the same as in Fig f4^). In the «-axis: 
1 s n/2, 2 eee n/ 3 , 3 ee r 2 /4, 4 ee r 3/4 , 5 ee r 3/5 and 6 ee r 4/5 



in Fig[4k) and the one from the ratios in Fig. |4j}). As expected 
the result for n < 2 is sensitive to the Model for the continuum 
which contributes for 40% to the moments. One can also note 
that using the moments from the data fit in [43] (Model 3), the 
result for n = 1 is: 



m c (m e )|i = 1289(8) MeV . 



(34) 



where the quoted error comes only from the change in a s given 
in Eq. (fT2l (some other sources of errors will be discussed 
later on). Though this result agrees with different determina- 
tions from A1i(0) 0371 138ll43l 14711 . one can note that its central 
value decreases when one increases the number of derivative n 
of the moments. The result only stabilizes versus the variation 
of n for n > 3 — 5, where an optimal result can be taken. For 
definiteness, we take n ^ 4 , where all Continuum Models give 
consistent predictions. Then, we deduce, from Fig|4^), in units 
of MeV: 



w c (m c )|o = 1263.7 



(1.3W3.5) B ,(4.9) a «(3.9) |l 
(4.4) C 2(4.7) C 3(3.1) C 4(1.7)„ ;) ,(35) 



where the central value is the average from different continuum 
models. It leads to the result in Table [2] 

The 1st error in Eq. (l35t is due to the different models for the 
continuum, the 2nd one to the value of a s given in Eq. (fl2l i. The 
3rd error is an estimate of higher order terms of PT assumed to 
be equal to the contribution of the a] one, while the 4th error is 
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an estimate of the effect of the subtraction point v by varying it 
from m c to M T and using the substitution (see e.g. Jj,Sl): 



a s (m c ) -> a s (v) X 1 1 - # — — log — 



(36) 



where P\ = -(1/2)(1 1 -2« y/3) for n /-flavours. The 5th and 6th 
errors are due respectively to the ones of the gluon condensates 
(a s G 2 ) and (g^fabcG^) estimated previously. The 7th error is 
due to the (G 4 ) condensates allowing it to move from fac=0.5 
(favoured value from our preceeding fit) to fac=l as defined in 
Eq. (1171 1, The last error is due to the experimental J/ip widths 
given in Table Q] 

We consider as a final value the one obtained from Alt(O) 
where both PT corrections are still small for the unit and di- 
mension 4 operators. Indeed, for the unit operator, the domi- 
nant correction is due to a s , which is about 14% for vVU(O) or 
for and which is about half of the one of A1i(0). Then, 
we may expect that the error induced by the organization of the 
PT series (fixed order, contour improved,...) discussed in [43] 
is smaller for vVU(O) though the PT series converges faster for 
vVf i(0) as on can notice in Eq. d32l . 

The result from the ratios of moments in Fig 03) is not very 
conclusive as the model-dependence of the result starts to dis- 
appear from the ratio of moments rys, but for these ratios the 
result increases with n. Then, we shall not retain the results 
from the ratios of moments for the charmonium channel. 

11. m c (m c ) from higher M„(Q 2 ) moments 

In the following, we shall extract Tfi c {m c ) from higher At n (4m^) 

and Ai„(8ml) moments. We show the results of the analysis 

respectively from the moments in Figs|5^) and|6^) and from the 

ratios of moments in Figs [5})) and|6j)). One can notice that in 

both cases the results from the moments present minimas versus 

n. 

The minimum is obtained from Ai\o(4m^) and from A1i5(8m^), 
which give in units of MeV: 



mc(mc)C,2 = 1261.9 



mMc)^ = 1260.9 



(0.7) ai (1.6) a 3(1.6) aP 4(0.4) /J 

(l.l) G :(1.0) G 3(1.7) G 4(3.0), t/) , 

(0.5) ar (1.6) a 3(1.6) af 4(0.4) // 
(1.0) G 2(0.7) G3 (1.3) G 4(2.6) eA7 ,, 

(37) 



which lead to the result in Table [2] The different sources of er- 
rors are the same as the ones discussed in Eq. (|35V The one 
from a\ here is due to the distance of the average of the a] con- 
tribution to the + assumed value of the coefficient in Eq. ( fTTT i. 
We have estimated the error due to the unknown a" (n > 4) to 
be equal to that of a? s . Eq. ( T37T i leads to the result in Table [2] 
One can also see in FigsO?) and[6t>) that the results from the ra- 
tios of moments increase with n. Though, the outputs obtained 
from the ratios of optimal moments are consistent with the ones 
from these moments and with the ones obtained in 111 where 
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Figure 5: a) Behaviour of m c {m c ) in units of MeV versus n from At„(4m7 ) 
moments and for the central value of (a s G 2 ) in Eq. 120) . We use fac=0.5 for 
the factorisation of (G 4 ) and Model 2 for the continuum. The colored region 
corresponds to the range of C3 values given in Eq. {TT}. b) Behaviour of m c (m c ) 
in units of MeV versus different ratios of moments r„/„+i(47n7 2 ) . The inputs 
and legends are the same as in Fig[3Ji). In the n-axis: 1 = 0/9, 2 = ("8/9, 3 = 

n/io, 4 = rg/io,. . . 



a judicious choice (small PT corrections) of these ratios have 
been used, we shall not consider these numbers in the final re- 
sults because of the absence of stabilities or minimas versus the 
variation of n. 



Table 2: Value of m c (m c ) from charmonium moments known to a] for Q 2 = and with 
an estimate of the cr s contribution for Q 1 ^ 0. 




Moments 


m c (m c ) [MeV] 


2 2 =0: 




M 4 


1263.7(10.3) 


Q 2 =4m 2 c : 




M w 


1261.9(4.5) 


<2 2 =8m 2 : 




M [5 


1260.9(4.0) 



7 The a 2 (resp a]) are relatively small i.e 7.4% (resp. 3.9%). 



12. Final value of m c (m c ) and Coulombic corrections 

Like in |ll], we approximately estimate the Coulombic correc- 
tion by working with the resummed expression of the spectral 
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Figure 6: a) Behaviour of m c (m c ) in units of MeV versus n from M„(8m^~) 
moments for the central value of (a s G 2 ) in Eq. {20) and of a s given in Eq. 112) . 
We use fac=0.5 for the factorisation of (G 4 ) and Model 2 for the continuum. 
The colored region corresponds to the range of cj, values given in Eq. 1111 . 
b) Behaviour of m c (m c ) versus different ratios of moments M n (Sm^ 2 ) in units 
of MeV. The inputs and legends are the same as in Fig|6ji). In the n axis: 
1 = 13/14, 2 = 1-13/15, 3 = 04/15, 4 = H4/16, 5 = n5/16>- ■ ■ 



function |48]0: 



'R-Acoui — -zv- 

2 1 — e 



(38) 



in an expansion in terms of x = Cpnajv where Cf = 4/3 
and v is the quark velocity defined in Eq. <j9j. This contribu- 
tion, which is of long-distance origin and proportional to the 
imaginary part (the wave function) of the two-point function, 
is induced by rescattering (Sommerfeld factor) of heavy quark 
pairs through the Coulomb potential above the cc threshold!! 
We add to this expression some PT QCD corrections. The 1st 
correction is the familiar (1 - ACfci s ) factor due to the quarko- 
nium annihilation through a single (transverse) virtual gluon. 
The 2nd t ype of corrections to order v and log v have been ob- 
tained in 1351 13611 up to order a 2 where the result is strictly ap- 
plicable near threshold Cfiia s < v <K 10 



8 However, we (a priori) expect that the non-relativisitc corrections will be 
relatively small here as we are working in the relativistic domain because m c is 
relatively light, while the final result corresponds to a relatively large Q 2 = 8m 2 
value. 

9 The Coulomb corrections arising from the bound states below the threshold 
can be safely neglected as the dispersion relation is applied above threshold 
(t > 4m 2 ) where the QCD expression of the spectral function from field theory 
(OPE) is used while its phenomenological expression is measured from the 
e + e~ — > hadrons data. 

'"However, according to Refs. [35, 36], these short-distance effects being 
specific for the single annihilation process involving QQ pairs are universal for 
|v| <K 1 regardless whether |v| is smaller or larger than Cpna s . 



We compare the value of the moments using the previous ex- 
pressions for the spectral function with the one obtained from 
PT theory including radiative corrections up to order a 2 . In the 
case Q 2 = %m 2 c and n = 15, where the most precise result is 
obtained, the corrections induced by the previous Coulombic 
contributions to the value of m c is about -1.3%E3E1 and gives: 

<Wc »; = ±16 MeV. (39) 

We consider this effect as another source of error rather than 
a definite shift on m c (jn c ) due to the fact that the role of the 
Coulombic effect in the sum rule analysis remains unclear 
1 3 igt] as the quark is still relativistic with a relatively large 
velocity: 



^(l + Q 2 /4m 2 Q ) In - 0.45 , 



(40) 



for large n — 15 and Q = %m Q . Indeed, this value of v would 
correspond to a momentum transfer between quark and anti- 
quark of about 1 GeV, where the effective potential differs from 
the Coulombic one [48] and where the sum rules are usually 
successfully applied. 

One can see from Table|2]that the estimate from different forms 
of the moments are consistent each other. We shall consider as 
a final estimate the most precise one from M l5 (%m 2 ), where 
the Coulombic correction obtained previously is also small. 
Adding this correction, we obtain: 



m c (m c ) = 1261(16) MeV , 
in excellent agreement with the one U] : 
m c (m c ) = 1261(18) MeV , 



(41) 



(42) 



obtained from a judicious choice of ratios of high moments hav- 
ing small PT and NP corrections. The previous results also im- 
prove earlier results obtained by the author to lower orders in 
this channel 15011 . 



13. Determination of mtimt) 

We extend the previous analysis to the bottomium systems, 
the following, we shall use the value: 



In 



a s (m b )\„ f=5 = 0.219(4) , 



(43) 



deduced from a s (m T ) in Eq. ( fl2l i. We shall use as experimen- 
tal inputs the T-family parameters in Table [3] using NWA and 
parametrize the spectral function above V? = (11.098 ± 0.079) 
GeV by its pQCD expression (Model 2), where the error in the 
continuum threshold is given by the total width of the T( 1 1 020) . 
We shall work with higher moments in order to minimize the 
contributions of the QCD continuum. We use moments known 
to order a] for Q 2 = 0, while for Q 2 ± 0, we have added the 
estimate of the a\ contribution given in Eq. (fTTT i. 



"Some further arguments justifying a much smaller value of these contri- 
butions can be found in 12411 - A much smaller effect of about 1 MeV is also 
obtained for the ratio of moments like has been found in [T|]. 

12 The effect on m c (m c ) from M l0 (4m 2 ) and A1 4 (0) are respectively 2 % and 
5%. 
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Table 3: Masses and electronic widths of the T family from PDG1O0. 
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Figure 7: a) Behaviour of mjfmj) in units of MeV versus n from M n (0) mo- 
ments for the central value of (a s G 2 ) in Eq. )20t and of a, given in Eq. 1121 . 
We use fac=0.5 for the factorisation of (G 4 > and Model 2 for the high-energy 
part of the spectral function, b) Behaviour of m;,(m/,) versus different ratios of 
moments A1„(0) in units of MeV. The inputs and legends are the same as in Fig 
[7^). In the rc-axis: 1 = rz/i, 2 = ra/4, 3 = 7-3/4, 4 = 7-3/5, 5 = r^/s,. . . 



Results from Ai„(0) 

We show the results from Ai„(Q) in Fig [7] where one can no- 
tice that the result is (almost) stable versus the variation of n 
for n =! 3 ~ 7 while for the ratios of moments, the stability is 
reached from /-4/g. At these values, the contribution of the QCD 
continuum is less than 29% of the total which is much less than 
the one for n = 1 where it is about 66%. This feature raises se- 
rious doubts on the accurate value of m b from this low moment 
M\{Q) given in the literature [37, 38, 47] due to the inaccuracy 
of the data in this high-energy region. From the moments, we 
obtain in units of MeV: 



m b {m b )t = 4169.6 



m b (m b )\ 4 /5 = 4166.3 



(1.9)^(4.1)^(2.7^ 

( 1 . 1 ) G 2 ( 1 .2) G 3 ( 1 .2) G 4( 10.6)^ . 

(4)^(7.1)^(16.6)^ 



4174 



Name 




Mass [MeV] 


rW,- [keV] 


Y(15) 




9460.30(26) 


1.340(18) 


4173 


Y(2S) 




lUUzj. 20(31 ) 


0.612(11) 


4172 


T(3S) 




10355.2(5) 


0.443(8) 


T(45) 




10579.4(1.2) 


0.272(29) 




T( 10860) 


10865(8) 


U. J l{ l ) 


4171 


Y(11020) 


1 1 ni q^s^i 


0.13(3) 
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Figure 8: a) Behaviour oirntirHb) in units of MeV versus n from M„(Amt") 
moments for the central value of (a s G 2 ) in Eq. HQ) and of a s given in Eq. 
1121 . We use fac=0.5 for the factorisation of (G 4 ) and Model 2 for the high- 
energy part of the spectral function. The colored region corresponds to the 
range of C3 values given in Eq. jilt , b) Behaviour of ~mi,(mt) versus different 
ratios of moments Ai„(4mb 2 ) in units of MeV. The inputs and legends are the 
same as in Fig[8ji). In the n axis: 1 = 0/8, 2 = r-j/g, 3 = rg/g, 4 a rs/io, . . . 



(0.6) C2 (0.4) C 3(0.4) C 4(1.9),, 



(44) 



giving the results in Table|4] One can notice that, at the optimal 
choice r4/5(0), PT corrections are large which induce larger PT 
errors than in the case of yVf4(0). The different sources of errors 
are similar to the case of charmonium. 

Results from A / („(4m^) 

The results from M„(4mf) are shown in Fig|8] where a stability 
versus the variation of n is obtained for n = 14, while for the ra- 
tios of moments, it is reached for rio/n(4m^) and rio/i2(4m^). In 
both cases, the errors due to the NP contributions and induced 
by the ± sign for the estimate of the a] coefficient and of the 
higher order a"- 4 are tiny (< 0.4 MeV) and can be neglected. 
We obtain in units of MeV: 



m b (m h )\ 



14 

4m? 



m b {m b )\2^ 



(45) 



4174.2(0.6)^(2.6)^(5.1)^ 
4178.6(4.2) ff ,(10.8y3.6), r/J , 
from which, we deduce the result in Table |4] 

Results from M„{%m 2 b ) 

The results from Ai„(8mfy and from the ratios of moments are 
shown in Fig [9] where stabilities versus the n-variations are re- 
spectively obtained for n — 9 ~ 11 and for ris/n, 06/17. Non 
perturbative corrections and the one due to the ± sign of the a 3 s 
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Table 4: Value of m\^n%) from bottomium moments known to a\ for Q 1 = and with an 
estimate of the a\ contribution for Q 1 + 0. 
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Figure 9: a) Behaviour of mj(mj,) in units of MeV versus n from M„{?Hn b ') 
moments for the central value of (a s G 2 ) in Eq. 1201 and and of a s given in Eq. 
1121 . We use fac=0.5 for the factorisation of (G 4 ). Model 2 for the high-energy 
part of the spectral function is used. The colored region corresponds to the 
range of values given in Eq. jilt , b) Behaviour of mj,(mj,) versus different 
ratios of moments M n (8mb~) in units of MeV. The inputs and legends are the 
same as in Fig[9^). In the n axis: 1 a ri3/14> 2 a H3/15. 4 a rn/15, . . . 



coefficient are also negligible (< 0.3 MeV). We obtain in units 
of MeV: 

mmt&i = 4175.1 (0.5)„ t (5.1) a 3(0.5)„ f 4 
(1.9)^(10)^, 
m b {m b )\Zi = 4170 - 9 (1.6)„,(9.5) a }(4.1) j* 

(7.2)^(3.6)^ . (46) 
Then, we deduce the result in Table [4] 



14. Final value of m b (m b ) and Coulombic corrections 

Here, we analyze the Coulombic corrections like in the case of 
charm. The ones for the moments are relatively large which 
are respectively 1.7%, 1.1% and 4% for M w (8m 2 b ), M l4 (4m 2 b ) 
and vV( 4 (0). The ones for the ratios of moments ?"i6/i7(8w|), 
'"io/n(4w?) and r4/5(0) are respectively 1.2, 2. 1 and 3.6 per 
mil, which are about one order of magnitude smaller. Among 
these different determinations the one from ri6/n(8m^) is the 
most precise. We consider it as our best final result: 

m b {m b ) = 4171(14) MeV. (47) 

It is informative to compare the previous result with the one in 
fl] (see Table from fl]) using some judicious choices of the 
ratios of moments having the smallest PT corrections and where 
the (G 4 ) contribution has not been included. Adding the errors 
±6 MeV due to the Coulombic, ±6 MeV due to the subtraction 



Moments 


m b (m b ) [MeV] 


Q 2 =0: 




> id 

M 


4169.6(11.3) 


n/5 


4166.4(18.6) 






M 14 


4174.2(5.8) 


T 10/11 


4178.6(12.1) 


2 =K : 




M i0 


4175.1(5.5) 


H6/17 


4170.9(12.7) 



point and ±4 MeV due to the a] contributions, the average from 
Table[5]becomes: 



TH b (m b ) = 4173(10) MeV, 



(48) 



which is in excellent agreement with the one obtained in Eq. 

63. 



Table 5: Corrected values of «7/,(m/,) from bottomiun moments known to 3-loops using 
some judicious choices of moments from Ref. 1 11]. The errors on my come respectively 
from the choice of the moments, a s , the data on the T family and the choice of the QCD 
continuum threshold. The ones due to the gluon condensates are negligible here. 



Mom 


m b (m b ) [MeV] 


2 =O: 




^2/3; r 2/A 


4160(4)(2)(3)(3) 


Q 2 =4ml : 




r S/9, r S/lO 


4177(2)03X3X6) 


Q 2 =%ml : 




r 13/14> f 13/15 


4183(2)(4)(2)(6) 


Average 


4173(4) 



15. Conclusions 

We summarize below the main results in this letter: 
We have explicitly studied in Section|5]the effect of the contin- 
uum model on the spectral function and found that this effect is 
large for Q 2 — low moments, which can only be evaded for 
moments Ai„>i-4(Q 2 )- This feature is naturally expected but 
raises the question on the errors induced by this model depen- 
dence in the determinations ofm c b from low-moments A / l„<2(0) 
used in the current literature. 

We have extracted the value of (a s G 2 ) in Section [7] and found 
the result in Eq. (1201 . This result confirms previous claims that 
the SVZ result underestimates the value of (a s G 2 ). We have 
not included in the analysis the most eventual short distance ef- 
fect of the D — 2 term advocated in 1151145311 which is dual to 
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the higher order terms of the PT series 115411 . However, like in 
different explicit analysis of some other light quark channels, 
the effect of this term might also be small and can improve the 
agreement between the QSSR predictions with the data or with 
some other determinations like lattice calculations. A future 
evaluation of this contribution is welcome but is beyond the 
scope of this paper. 

We have re-extracted the gluon condensate (g 3 f a bcG 3 ) and ob- 
tained its value in terms of the instanton radius p c in Eq. ( f25l l. 
This value agrees within the error with the one in QJJ] but is 
smaller than the estimate from the DIG approximation p c - 5 
GeV 

During the determinations of these condensates, our analysis 
prefers the value fac=0.5 of the D = 8 (G 4 ) condensates which 
supports the modified factorisation proposed in ll45ll using a 
1 IN approach. 

We have re-estimated the MS running masses ~m c j, to order a\ 
and including the (G 4 ) condensate contributions in the OPE. 
Optimal results from different moments lead to the final values 
in Eqs. ( HTT i and (PTTT i. These results confirm the recent results in 
IH obtained from judicious choices of ratios of moments with 
small PT corrections and where the contributions of the D = 8 
(G 4 ) condensates have not been included. They also improve 
older results in [50] obtained at lower orders with larger errors. 
These results are also comparable with the ones in the existing 
literature using different methods iH |37l l43l WA |55l l56ll . 
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